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Disorder induced brittle to quasi-brittle transition in fiber bundles 
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PACS 64.60.av - General studies of phase transitions: Cracks, sandpiles, avalanches, and earth- 
quakes 

PACS 46.50.+a - Fracture mechanics, fatigue and cracks 

Abstract. - We investigate the fracture process of a bundle of fibers with random Young modulus 
and a constant breaking strength. For two component systems we show that the strength of 
the mixture is always lower than the strength of the individual components. For continuously 
distributed Young modulus the tail of the distribution proved to play a decisive role since fibers 
break in the decreasing order of their stiffness. Using power law distributed stiffness values we 
demonstrate that the system exhibits a disorder induced brittle to quasi-brittle transition which 
occurs analogously to continuous phase transitions. Based on computer simulations we determine 
the critical exponents of the transition and construct the phase diagram of the system. 



I . Introduction. — The fracture of heterogeneous ma- 
' O teriais is a very important scientific problem with a broad 
G spectrum of technological applications ■ During the 
Qpast decade large efforts have been devoted to obtain a 
!— deeper understanding of the role of disorder in fracture 
^ processes. On the one hand, the problem is very important 
from a technological point of view in order to be able to de- 
00 sign novel type of composite materials with high strength 
■^j"at a low weight. On the other hand, due to the decisive 
C role of disorder, fracture phenomena address interesting 

challenges also for statistical physics [TJ[2] ■ 
^^O ' One of the most important modeling approaches to the 
< ^ > fracture of heterogeneous materials is the fiber bundle 
^^rhodel (FBM) [2H3- In the framework of the model the 
• • sample is discretized in terms of a parallel bundle of fibers. 
m The system is typically loaded parallel to the fiber direc- 
tion by slowly increasing the external load or the deforma- 
$— <tion. Assigning appropriate physical properties to fibers 



such as Young modulus, breaking strength, and Theolog- 
ical behaviour, several important features of the fracture 
process can be reproduced. During the last two decades 
FBMs have provided a very useful insight into the fracture 
of heterogeneous media both on the macroscopic and mi- 
croscopic length scales [2H§]. A very important common 
feature of all these investigations is that the Young modu- 
lus of fibers is assumed to be constant, and the disorder of 
the material is captured solely by assigning random break- 
ing thresholds to the fibers. However, it is well known that 
composite materials made up as a mixture of several in- 
gredients, can also have a higher degree of disorder where 



not only the strength but even the Young modulus of local 
material elements can vary [2]. 

The goal of the present Letter is to capture the het- 
erogeneity of the stiffness of elements and investigate the 
emerging complex fracture process in the framework of 
fiber bundle models. We focus on aspects relevant from 
the viewpoint of statistical physics and show that the sys- 
tem has several interesting novel features. Assuming a 
constant strength for fibers, under global load sharing the 
strain is everywhere the same in the bundle, hence, the 
fibers break in the decreasing order of their stiffness val- 
ues. It has the consequence that the tail of the disor- 
der distribution plays a crucial role in the breaking pro- 
cess. For two component systems analytical calculations 
revealed that the macroscopic strength of the mixture is 
always lower than the strength of its components. As the 
main outcome, our investigations showed that in a bun- 
dle of power law distributed stiffnesses a disorder induced 
transition occurs from a perfectly brittle phase where the 
first fiber breaking triggers the immediate collapse of the 
bundle to a quasi-brittle phase where precursors emerge 
before failure. Based on analytical calculations and com- 
puter simulations we construct the phase diagram of the 
system and determine the critical exponents. 

FBM with Random Young modulus. fn the 

classical FBM a bundle of N fibers is considered which all 
have a linearly elastic behavior characterized by the same 
Young modulus E. The fibers can sustain a finite load, 
i.e. when the load on them exceeds a threshold value o~th, 
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the fibers break irreversibly. To capture the heterogeneity 
of the local physical properties of materials it is assumed 
that the strength of fibers a t h is a random variable with 
a probability density p{a t h) and a cumulative distribu- 
tion P(ath) = jn th p[x)dx. After a fiber breaks, its load 
has to be overtaken by the remaining intact ones which 
introduces interaction between the fibers. For the load re- 
distribution two limiting cases are usually analyzed: equal 
load sharing (ELS) means that all the intact fibers share 
the same load irrespective of their distance from the failed 
one [cHHUIi] . ELS implies that there is no stress concentra- 
tion in the bundle. In the opposite limit of localized load 
sharing, the load is redistributed over the close vicinity 
of broken fibers leading to strong overloads around failed 
regions [6]. 

Assuming equal load sharing the macroscopic constitu- 
tive equation a(s) of the system can easily be obtained 
analytically. Loading the bundle parallel to the fibers' 
direction, the strain s is everywhere the same in the sys- 
tem. At a given e all the intact fibers keep the load Ee 
and their fraction is given by [1 — P(Ee)}, since the fibers 
with breaking thresholds ath < Ee have already failed. It 
follows that the constitutive equation takes the form 



cr(e) =Ee[l- P{Ee)\ 



(1) 



It has been shown in the literature that for a broad class 
of disorder distributions p(a t h) the constitutive curve cr(e) 
has a quadratic maximum whose position and value define 
the critical deformation e c and critical stress a c of the 
system, respectively PHSHZHU]- 

In the present paper we modify the classical FBM by 
introducing randomness for the Young modulus of fibers 
described by the probability density function p(E) over 



the interval E v 



<E<Er, 



For simplicity we assume 



that the strength of fibers is constant ath = const so that 
the random Young modulus is the only source of disorder. 

Constitutive equation. In order to derive the constitu- 
tive equation we assume that the system is loaded through 
rigid bars which ensure that the deformation e is every- 
where the same in the bundle. The random Young modu- 
lus of fibers has the consequence that at a given value of e 
the fibers carry different loads at — E%e, where Ei denotes 
the Young modulus of fiber i (i — 1, . . . ,N). Slowly in- 
creasing the externally imposed deformation e those fibers 
will first reach the constant breaking threshold ath, which 
have the highest Young modulus, i.e. for which a t h — E^e 
hold. It follows from the above arguments that in the 
bundle fibers break in the decreasing order of their Young 
modulus and at a given strain e those fibers are broken 
whose Young modulus exceeds ath/z- 

The load kept by the fibers having Young modulus be- 
tween E and E + dE reads as Eep(E)dE which has to be 
summed up from E m i n to the upper limit of intact fibers 
ath/s- Hence, the the constitutive equation of a fiber bun- 
dle with random Young modulus and a constant breaking 



threshold can be cast into the generic form 
a(e) = £ / Ep{E)dE. 

J E min 



(2) 



To understand the mechanical behavior of the bundle it is 
instructive to consider some limiting cases of the loading 
process: at small deformations e — > 0, the upper limit of 
integration in Eq. (0) goes to infinity so that we have 



a(e ->■ 0) w e 



p{E)EdE = (E) e, 



(3) 



E min 



where (E) denotes the average value of the Young modu- 
lus of fibers. The result implies that the system displays 
linear behavior in the limit of small strains, where the 
macroscopic Young modulus of the bundle is equal to the 
average Young modulus of the fibers. At high deforma- 
tions e — y oo the upper limit of integration goes to zero so 
that a(e — > oo) = indicating the breaking of all fibers. 
Between the two limits, the constitutive curve can have a 
maximum whose position e c is obtained from the equation 



da 
ds 



(rth/ec 



,,[E)EdE-[ °-f- ) plrra,/:,.) = (). (4) 



defining the critical strain of the bundle. Under stress 
controlled loading the bundle fails catastrophically when 
reaching the maximum stress a c = a(e c ), which is the 
macroscopic strength of the bundle. 

Mixture of two components. Recently, it has 
been shown that forming a composite of hard and soft 
components can increase the strength of the material. 
One example of such two-component composites where 
strengthening occurs is a so-called double-network (DN) 
gel, where a network of brittle polyelectrolyte gel and flex- 
ible polymer chains are mixed [lOllllj . As an application 
of our modeling approach, first we consider a simple two- 
component mixture composed of two subsets of fibers with 
different Young modulus E\ and £2 where E\ < E^. The 
fractions of the two subsets arepi andp2, where p\-\-p2 = 1 
holds. The main parameters of the model are r = E2/E1, 
and k = pijp\ with ranges 1 < r < 00 and < k < 00, 
respectively. 

Starting from the generic expression Eq. ([2]) it follows 
that before the first breaking event the bundle has a lin- 
early elastic behavior with the constitutive relation 



a = s(EiPx + E2P2) = eExp^l + rk), 



(5) 



where the average Young modulus has the expression 
(E) = p\E\ + p2E 2 . Increasing the external load, first 
fibers of the subset with the higher Young modulus £2 
break. This breaking occurs instantaneously, i.e. all the 
fibers of the subset are removed at once when a t h — E 2 e 
is reached. For the stress at the first breaking a\ we can 
write (Ji = E X (E) = a th (pi/r +p 2 )- 
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Fig. 1: (Color online) Macroscopic response of the two- 
component system for two different values of r. The solid 
red line represents the constitutive curve under stress con- 
trolled loading while the blue dashed line shows the behavior 
under strain controlled loading. The value of the parameters 
Ei = l,pi = 0.6, p2 = 0.4, k = 0.66, a t h = 10 are fixed. In (a) 
for r = 5 the stability condition Eq. ([B]) holds so the system 
stabilizes after the first collapse and the macroscopic strength 
of the bundle is determined by the second subset. In (b) for 
r = 2 the stability condition fails so that macroscopic failure 
occurs right at the breaking of the first subset. 



Performing strain controlled loading, after the first col- 
lapse (of the fibers with the higher Young modulus E2) 
the strain of the system e remains the same while the 
stress drops down as a result of the breaking by an amount 
Act = <JthP2 to the new value a[ = (T t hPi/r. After the 
stress drop, we can keep loading the system until we 
reach the second collapse, where all fibers with the lower 
Young modulus E\ break. This event occurs at the strain 
£2 = &th/Ei and stress a 2 = £2 (E') = uthPi values. Com- 
paring a± and 02 we can check whether the maximum load 
capacity a c of the system is determined by the collapse of 
the first or the second subset. Based on the above results 
it can be derived that if the parameters of the model r 
and k fulfill the condition 



k < 1 - l/i 



(6) 



the critical load <r c of the system is determined by the 
second collapse a c = cr 2 , otherwise, by the collapse of the 
first subset with the higher Young modulus a c = o\. 

Under stress controlled loading conditions, the sys- 
tem follows the same constitutive curve a(e) as under 
strain controlled loading until the first collapse. Assum- 
ing a c — (T2, at this point the stress is kept constant and 
the bundle gets elongated by an amount Ae = kath/Ei 
to the new value e' — othjEi + ka t h/Ex where the re- 
maining fibers will be able to sustain the load. Figure [1] 
presents the macroscopic response of the two-component 
system for two different parameter sets: in (a) the condi- 
tion k < 1 — 1/ r holds so that stability is retained after the 
failure of the stiffer subset of fibers, while in (b) macro- 
scopic failure occurs right after the first breaking. We 
indicated the path of both the strain and stress controlled 
loading. 

The most important characteristics of the system is the 
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Fig. 2: (Color online) Strength of the two-component system 
normalized by the breaking threshold a t h as a function of (a) 
the ratio r = E2/E1 with fixed k = 2/3 and (b) the ratio 
k = P2/P1 with fixed r = 2. From both (a) and (6) it can be 
seen that the single component system is stronger than the two 
component one. 



macroscopic strength a c . Based on the above results we 
can express the ratio of the fracture strength of the bundle 
er c and of the breaking threshold ath as function of the 
composition parameters k and r. For fixed values of k we 
obtain 



<J c /vth{r) = 



P2, 



Pi, 



1 < r < 

l-k ^ 



1 . 

l-k- 



and for fixed r the calculations yield 



c c /W(fc) 



1/(1 + *), 

l+rk 
r-\-rk 5 



< k < 1 



1 



< k. 



(7) 



(8) 



The results are plotted in Fig. [2j The main outcome of the 
calculations is that in the framework of the fiber bundle 
model with infinite range of interaction, i.e. global but 
not equal load sharing ensured by the rigid loading bars, 
the single component system has always higher strength 
a c than a two-component mixture. The worst situation 
with the lowest strength is obtained for the parameter 
combination k = 1 — 1/r. The results are in agreement 
with Ref. [TU] where it was shown that the strengthening 
in two component mixtures only occurs if a large enough 
crack exists in the sample before the loading starts. 

Disorder induced brittle to quasi-brittle transi- 
tion. — As the next step of the investigations we consider 
the case when the Young modulus of fibers has a continu- 
ous distribution p{E) over the interval E m i n < E < E max . 
Starting from the generic form of the constitutive equation 
Eq. ([5]) it can be shown that for a uniform distribution of 
E the breaking of the first fiber with the highest Young 
modulus results in immediate catastrophic collapse of the 
entire bundle. The reason is that for the uniform dis- 
tribution a relatively large fraction of fibers have Young 
modulus in the vicinity of the largest value E max , and 
hence, the load increment created by the first breaking 
can trigger a catastrophic avalanche. For the Weibull dis- 
tribution such catastrophic collapse does not occur due to 
the exponential tail of the distribution. These arguments 
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Fig. 3: (Color online) (a) Macroscopic behavior of the bun- 
dle with power law distributed Young modulus for different 
values of the exponent a and fixed values of E, n i n — 0.5 and 
E max — 1. (b) Phase diagram of the fiber bundle with power 
law distributed Young modulus on the r-a plane. The phase 
boundary between the perfectly brittle and quasi-brittle macro- 
scopic behavior is formed by the curve of r c (a). 



show the importance of the shape of the stiffness disorder 
in the vicinity of the upper bound, hence, in the follow- 
ing we consider a power law distribution p{E) = CE~ a 
for the Young modulus, where the functional form can be 
controlled by varying the exponent < a. In the limiting 
case a = we recover the uniform distribution, while in- 
creasing a the function p(E) decreases faster making the 
majority of fibers less stiff. 

The constitutive equation a(e) of the system can be 
determined by substituting the distribution p(E) after 
proper normalization into the generic form of Eq. (|2"|) 



CTa(e) 
(71 (s) 

£72 (e) 



1 - a 



1 



* " rjrnax — & rnin 

~ i&th E m i n e) . 

In ■ 



\?th e 



E min £ ) > 



(9) 



E min 



in 



e — In (e) e 



where cr a (e), a\(e), (72(e) correspond to a generic a differ- 
ent from I and 2, to the special case of a = 1 and a = 2, 
respectively. Figure \3^a) shows representative examples 
of the constitutive curve a(e) for different values of the 
exponent a with fixed E m i n = 0.5, E max = 1. It can 
be observed that at low values of a the constitutive curve 
(7(e) has a sharp maximum, i.e. 17(e) has a linear behavior 
up to the maximum followed by immediate collapse which 
indicates a perfectly brittle behavior. However, at high a 
values (7(e) develops a quadratic maximum, i.e. nonlinear 
behavior occurs before the maximum is reached. It follows 
that in this parameter range, where the stiffness distribu- 
tion decreases rapidly, the system becomes quasi-brittle 
being able to suffer several avalanches of fiber breaks be- 
fore collapsing in a catastrophic avalanche. In this case the 
macroscopic failure of the system is preceded by precur- 
sors which provide important signals of the imminent fail- 
ure event. The boundary of the two phases can be found 
by analyzing the derivative of the constitutive curve da j de 




Fig. 4: (Color online) The burst size distribution for power law 
decreasing stiffness disorder for 10 6 fibers averaged over 1000 
samples with different a values and a fixed r = 0.5. Approach- 
ing the phase boundary a —¥ a c — 3 a crossover occurs from 
the exponent 5/2 to a lower one 3/2. 



at the point of the first fiber breaking e = <Jth/ E max . Dif- 
ferentiating Eq. (TPTj) we get 



d<T a 
de 

da 1 

de 

da 2 

de 



E n 



E n 



E n 



1 - a 

'2 - a 

r 

- 1 ' 
r 

'r - 1 



a — 1 — r 2 
1 -r 1 -' 



[l + lnr] , 



. (10) 

(11) 
(12) 



where the parameter r — E m i n /E max can vary in the 
range < r < 1. If the derivative is negative the bun- 
dle is perfectly brittle and collapses without precursors, if 
however it is positive da / de\ C!th / E max > precursors can 
be observed and a quasi-brittle behavior emerges (see also 
Fig- El)- Analytic calculations show that the stability fea- 
tures are determined by the range of Young modulus r 
and by the exponent a: For the interval < a < 1 the 
system collapses at the instant of the first fiber breaking 
regardless of the value of r. This also means that for uni- 
formly distributed Young modulus a — no stability can 
be obtained. For a > 1 the stable regime of r is r > r c 
where the critical value r c as a function of a can be cast 
in the form 

r c (a) = (a-l) 1 /(2-). (13) 

The result implies that at a given value of the exponent a 
the range of Young modulus values has to be small enough 
to stabilize the system, or if we fix the range r the distri- 
bution has to decay fast enough to obtain stability. The 
results are summarized in Fig.[3j6) which provides a phase 
diagram of the system on the r-a plane. The boundary 
between the phases of perfectly brittle and quasi-brittle 
behaviors is provided by the function r c (a) Eq. (|13[) . 

Microscopic dynamics - statistics of avalanches. Un- 
der stress controlled loading, after each breaking event 
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Fig. 5: (Color online) (a) Average avalanche size as a func- 
tion of a for several different system sizes N. As -/V increases 
the position of the maximum a c (N) shifts towards the critical 
point of the infinite system a c (oo) determined analytically, (b) 
Finite size scaling of the average cluster size. The good quality 
collapse was obtained with the exponents 7 = 2 and v — 3. 



the load of the broken fiber gets redistributed over the in- 
tact ones which can induce additional breakings and finally 
can even trigger an entire avalanche of failing fibers. The 
size of the avalanche A is defined as the number of fibers 
breaking in the avalanche. The emergence of avalanches 
is the direct consequence of the quasi-brittle macroscopic 
response in FBMs which can also be used to forecast the 
imminent failure event. It is important to emphasize that 
in our FBM where the fibers have random Young modu- 
lus, the load carried by the fibers is different even if we 
assume global load sharing. 

In order to study the microscopic avalanche dynamics 
computer simulations were carried out at a fixed range of 
Young modulus with r = 0.5 which yields a critical value 
of a c = 3 for the infinite system from Eq. (|T3"]) . To char- 
acterize the statistics of avalanches we determined their 
size distribution p(A) for several values of the exponent a 
inside the stable regime between the limits a c < a < 10 
as it is represented by the arrow in Fig. It can be 

observed in Fig. 0] that the burst size distribution has a 
power law behavior 



P (A) ~ A"- 



(14) 



for all the parameter values considered. It has to be 
emphasized that as we approach the phase boundary a c 
from the quasi-brittle phase, the distribution p( A) shows a 
crossover from a higher to a lower exponent: Far from the 
phase boundary the value r = 5/2 is obtained in agree- 
ment with the usual exponent of FBMs with equal load 
sharing [2J[U[!5]. As the response becomes more brittle 
a a c {r) the exponent switches to a lower value r = 3/2. 
This crossover is similar to what is observed in classical 
FBMs with solely strength disorder when increasing the 
lower bound of the range of strength values towards the 
critical one of immediate collapse [5,6 . However, in our 
case the range of Young modulus r is fixed, while the func- 
tional form of the disorder distribution is changed. The 
results indicate that the higher degree of brittlcness im- 
plies a lower exponent r of the burst distribution due to 



the higher frequency of avalanches of larger size. 

In order to reveal the nature of the transition from the 
quasi-brittle phase to the one of perfectly brittle behavior 
with catastrophic collapse, we extended the simulations 
to a values inside the brittle phase. We determined the 
average size of avalanches (A) as a function of a for several 
different system sizes N. The average size of bursts (A) is 
calculated as the ratio of the second and first moments of 
the burst size distribution skipping the largest avalanche 



<A) = 



(15) 



It can be observed in Fig. (5][a) that the average avalanche 
size has a maximum as a is varied. Increasing the system 
size N, the maximum gets more and more peaked and its 
position shifts to lower a values. Left from the maximum 
(A) decreases rapidly and goes to zero marking the per- 
fectly brittle response in this a regime. It follows that 
the position of the maximum of (A) can be identified as 
the pseudo-critical point a c (N) of the finite size system 
which converges towards the critical point of the infinite 
system a c (oo) determined analytically in the previous sec- 
tion. Assuming that the system has a continuous phase 
transition from the quasi-brittle to the brittle phase when 
a decreases, the finite size scaling form 



a c (N) = a c {oo)+BN- 1/u 



(16) 



has to hold, where v is the correlation length exponent of 
the transition [12] . Figure [SJ&) presents that by rescaling 
the curves of (A) (a, N) by an appropriate power of the 
system size N taking into account the functional form Eq. 
(|16l) of a c (N), the curves of Fig. [5ja) can be collapsed on 
top of each other. The high quality collapse implies that 
(A) (a, N) has the scaling structure 



(A) (a, N) = N^/"F^((a ~ a c )N l ' v ), 



(17) 



where 7 is the susceptibility exponent of the brittle to 
quasi-brittle transition and denotes the scaling func- 
tion [12] . Based on the data collapse analysis the values 
7 = 2 and v = 3 were obtained. 

Perfectly brittle behavior of the system implies that no 
damage can be accumulated prior to failure. However, in 
the quasi-brittle phase a finite Nb number of fibers break in 
avalanches before the catastrophic avalanche. Hence, we 
introduce the order parameter of the system as the frac- 
tion of fibers breaking prior to failure pb — Nb/N which is 
zero in the brittle and has a finite value in the quasi-brittle 
phase. For the infinite system pb can be obtained analyt- 
ically as pb = 1 — P{<Jth/£c), where P is the cumulative 
distribution of the Young modulus and e c is the critical 
deformation of the system obtained from Eq. (f5]) . Finally, 
Pb can be cast into the form 



Pb 



( a _ l)(a-D/(2-«) 



(18) 
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(a - a c ) ■ N l < v 

Fig. 6: (Color online) Inset: fraction of fibers pi, breaking before 
final failure as a function of a for different system sizes together 
with the analytic solution Eq. (|18fl . Main panel: Finite size 
scaling of pt- The good quality collapse was obtained with 
f3 — 1 and v = 3. 



which has a power law dependence on the distance from 
the critical point p^ ~ (fl — a c (f))^. The value of the order 
parameter exponent is = 1. The inset of Fig. H2 presents 
Pb as a function of a obtained by simulations for differ- 
ent system sizes together with the analytic solution Eq. 
(fF8l . The main panel of the figure demonstrates the good 
quality collapse of Pb(a, N) assuming the scaling structure 

Pb (a, N) = N~^ u F (2 \(a - a c )N^ u ). (19) 

The best collapse was obtained with the exponents (3 = 1 
and v = 3 showing the consistency of the results. The 
scaling function is denoted by . It is interesting to note 
that the critical exponents /3, 7, and v fulfill the hyper- 
scaling relation /3 + 7 = v. 

Summary. — We investigated the fracture process of 
heterogeneous materials in the framework of a fiber bun- 
dle model. As a novel feature of the model we assumed 
that the source of disorder is the randomness of the Young 
modulus of fibers, while their breaking strength is kept 
constant. We showed by analytical and numerical calcula- 
tions that the system exhibits several interesting features 
both on the micro and macro scales. For two component 
systems we demonstrated that the strength of the mixture 
is always lower than that of the individual components. 

To obtain a deeper understanding of the breaking pro- 
cess, a power law distribution was introduced for the 
Young modulus defined over a finite range. We showed 
analytically that varying the range of the Young modulus 
and the value of the exponent of the disorder distribution 
the system exhibits a transition from a perfectly brittle 
phase where the breaking of the first fiber triggers the im- 
mediate collapse of the system to a quasi-brittle one where 



macroscopic failure is preceded by avalanches of breaking 
events and by a non-linear macroscopic response. Based 
on analytic calculations and computer simulations we con- 
structed the phase diagram of the system and we showed 
that the brittle to quasi-brittle transition occurs analo- 
gously to continuous phase transitions. The critical expo- 
nents of the transition, determined by finite size scaling 
analysis, fulfill the hyper-scaling relation f3 + 7 = v. It 
is a crucial feature of the system that due to the global 
load sharing there is no spatial correlation between fiber 
breakings, i.e. broken fibers nucleate completely randomly 
all over the bundle. This is the reasons why the brit- 
tle to quasi-brittle transition shows some similar features 
to percolation. However, it has to be emphasized that 
the macroscopic failure of the bundle has nothing to do 
with the appearance of a spanning cluster of broken fibers 
which is also indicated by the fact that the fraction of bro- 
ken fibers pb at global failure is relatively small. For future 
investigations, it would be very interesting to consider the 
competition of two sources of disorder, i.e. to investigate 
the case where both the Young modulus and the breaking 
strength of fibers are random. Work in this direction is in 
progress. 
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